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1. Introduction 

The puzzle of the hypothetical dark matter 
(DM) in the Universe still remains unresolved . 
While there may well be more than one type of 
DM, arguments from large structures suggest that 
a large, and presumably dominant, fraction of 
DM in the Universe is made of massive particles 
which at the time of entering the epoch of mat- 
ter dominance would be already non-relativistic, 
or cold. From the particle physics point of view, 
cold DM (CDM) could most plausibly be made 
of so-called weakly interacting massive particles 
(WIMPs). 

There exist several interesting WIMP candi- 
dates for CDM that are well-motivated by the 
underlying particle physics. The neutralino is 
considered by many a "front-runner" by be- 
ing perhaps the most natural WIMP: it comes 
as an unavoidable prediction of supersymmetry 
(SUSY) . The axion is another well-motivated can- 
didate || . But by no means should one forget 
about some other contenders. While some old 
picks (sneutrinos and neutrinos with mass in the 
GeV range) have now been ruled out as cosmolog- 
ically relevant CDM by LEP, axinos (SUSY part- 
ners of axions) have recently been revamped ||]. 
Another well-known candidate is the fermionic 
partner of the graviton, the gravitino. Most of 
the review will be devoted to the neutralino but I 
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will also comment about recent results regarding 
axinos. 

Neutrinos, the only WIMPs that are actually 
known to exist, do not seem particularly attrac- 
tive as DM candidates. It has long been believed 
that their mass is probably very tiny, as suggested 
by favored solutions to the solar and atmospheric 
neutrino problems, which would make them hot, 
rather than cold DM. This picture has recently 
been given strong support by first direct evidence 
from Superkamiokande for neutrinos' mass 
While the new data only gives the [i — t neu- 
trino mass-square difference of2.2xl0~ 3 eV 2 , it 
it very unlikely that there would exist two mas- 
sive neutrinos with cosmologically relevant mass 
of 5 to 40 eV and such a tiny mass difference. 

The strength with which the WIMP candidates 
listed above interact with ordinary matter spans 
many orders of magnitude. For the neutrali- 
nos it is a fraction (~ 1CU 2 — ~ 10~ 5 ) of the 
weak strength. Interactions of axions and axi- 
nos are suppressed by (mw/ fa) 2 ~ 10~ 16 , where 
f a ~ 10 10 GeV is the scale at which the global 
Peccei-Quinn U(l) symmetry is broken. Interac- 
tions of gravitinos and other relics with only grav- 
itational interactions are typically suppressed by 
(mw/Mpi anc k) 2 ~ 10 -33 . One may wonder how 
such vastly different strengths can all give the 
relic abundance of the expected order of the crit- 
ical density. The answer lies in the different ways 
they could be produced in the early Universe: the 
neutralinos are mainly produced "thermally" : at 
some temperature they decouple from the thcr- 
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mal bath. But in addition there exist also several 
mechanisms of "non-thermal" production. This 
is how CDM axinos, gravitinos can be produced 
(in addition to "thermal" production), as well as 
Wimpzillas ||. 

It is obvious that WIMPs do not necessarily 
have to interact only via weak interactions per 
se. WIMPs are generally required to be electri- 
cally neutral because of stringent observational 
constraints on the abundance of stable charged 
relics. On the other hand, they could in princi- 
ple carry color charges. (For example, a stable 
gluino above 130 — 150 GeV ||, or in an experi- 
mentally allowed window 25 — 35 GeV [Q, could 
still be the lightest SUSY particle (LSP) although 
its relic abundance would be very small ||.) I n a 
halo they would exist as neutral states by form- 
ing composites with gluon or quarks. We can see 
that generally WIMPs are expected to have sup- 
pressed effective couplings with ordinary matter. 
Otherwise, they would dissipate their kinetic en- 
ergy. 

What is the WIMP relic abundance Q x h 2 — 
Px/Pcrit Hi? Current estimates of the lower 
bound on the age of the Universe lead to 
^total/i 2 < 0.25. Recent results from rich clus- 
ters of galaxies and high-redshift supernovae type 
la imply fimatter — 0.3. The Hubble parameter is 
now constrained to 0.65 ± 0.1. Big Bang nucle- 
osynthesis requires fibaryon^ 2 < 0.015. Assuming 
that most matter in the Universe is made of CDM 
WIMPs, one therefore obtains 

0.1 <n x h 2 < 0.15 (1) 

or so. At the very least, requiring that a dominant 
fraction of CDM is located only in galactic halos 
gives Vt x h 2 > 0.025. 

2. Neutralinos 

The DM candidate that has attracted perhaps 
the most wide-spread attention from both the 
theoretical and experimental communities is the 
neutralino. It is a neutral Majorana particle, the 
lightest of the mass eigenstates of the fermionic 
partners of the gauge and Higgs bosons: the bino, 
wino and the two neutral higgsinos. If it is the 
lightest SUSY particle, it is massive and stable, 



due to assumed R-parity. A perfect candidate for 
a WIMP! 

The neutralino is a well-motivated candidate. 
It is an inherent element of any phenomenologi- 
cally relevant SUSY model. Being neutral, it is 
a natural candidate for the LSP (although one 
should remember that this is often only an as- 
sumption); it couples to ordinary matter with a 
weak-interaction strength (reduced by mixing an- 
gles) which is within the range of sensitivity of 
present-day high energy, as well as dark matter, 
detectors. Finally, as a bonus, it naturally gives 
Vl x h 2 ~ 1 for broad ranges of masses of SUSY 
particles below a few TeV and for natural ranges 
of other SUSY parameters. 

Much literature has been devoted to the neu- 
tralino as DM, including a number of comprehen- 
sive and topical reviews (see, e.g., Refs. ||,|l(|). 
Here I will only summarize the main results and 
comment on some recent developments and up- 
dates. 

Neutralino properties as DM and ensuing im- 
plications for SUSY spectra are quite model de- 
pendent but certain general conclusions can be 
drawn. First, its cosmological properties are very 
different depending on a neutralino type. The 
relic abundance of gaugino-like (mostly bino-like) 
neutralinos is primarily determined by the (light- 
est) sfermion exchange in the LSP annihilation 
into //: ilh 2 oc m~/m x . In order to have fi x ~ 1, 

the lightest sfermion cannot be too light (below 
~ 100 GeV) nor too heavy (heavier than a few 
hundred GeV) jllj which is a perfectly natural 
range of values. 

On the other hand, higgsino-like V s are 
strongly disfavored. Firstly, due to GUT relations 
among gaugino masses, higgsinos correspond to 
rather large gluino mass values, above 1 TeV, 
which may be considered as unnatural [|ll) . Fur- 
thermore, higgsino-like neutralinos have been 
shown to provide very little relic abundance. For 
m x > mz,mw,mt the x pair- annihilation into 
those respective final states (ZZ, WW, tt) is 
very strong. But both below and above those 
thresholds, there are additional co-annihilation 
|l2] ] processes of the LSP with xt and which 
are in this case almost mass-degenerate with the 



4 



LSP. Co-annihilation typically reduces Vt x h 2 be- 
low any interesting level although it has been re- 
cently argued that the effect of co-annihilation is 
not always as strong as previously claimed |13| . 
Higgsino-like LSPs are thus rather unattractive, 
although still possible DM candidates, especially 
in the large mass (m x > 500 GeV) regime. One 
also arrives at the same conclusions in the case 
of 'mixed' neutralinos composed of comparable 
fractions of gauginos and higgsinos. This is be- 
cause, even without co-annihilation, in this case 
the neutralino pair-annihilation is less suppressed 
and one invariably finds very small fl x there p"l) . 

Remarkably, just such a cosmologically pre- 
ferred gaugino type of neutralino typically 
emerges in a grand-unified scenario with addi- 
tional assumptions that the mass parameters of 
all spin-zero particles are equal at the unification 
scale ~ 10 16 GeV. What one finds there is that 
the lightest bino-like neutralino emerges as essen- 
tially the only choice for a neutral LSP Jl4|,[l5| . (It 
is still possible to find cases with a higgsino-like 
LSPs but they are relatively rare.) Furthermore, 
in order for Q, x h 2 not to exceed one or so, other 
(most notably sfermion) SUSY partner masses 
have to be typically less than 1 TeV 
Thus our phenomenological expectations for low- 
energy SUSY to be realized roughly below 1 TeV 
are nicely consistent with a bino-like neutralino 
as a dominant component of DM in the Universe. 

What about the neutralino mass? In princi- 
ple, it could be considered as a nearly free pa- 
rameter. It would be constrained from below 
by the requirement Q, x h 2 < 1 to lie above some 
2 — 3 GeV jl6| which is a neutralino version of the 
Lee- Weinberg bound [Q. This is because, in this 
mass range, the neutralino may be viewed as a 
massive neutrino with somewhat suppressed cou- 
plings. Much stronger bounds are in many cases 
provided by LEP. Unfortunately, collider bounds 
are rather model dependent and are no longer 
provided for the general SUSY model normally 
considered in DM studies. Nevertheless, one can 
reasonably expect that m x is now ruled out below 
roughly 30 GeV, except in cases when the sneu- 
trino is nearly degenerate with the LSP in which 
case the lower bound may disappear altogether. 

A rough upper bound of 1 TeV follows from 



the theoretical criterion of naturalness: expecting 
that no SUSY masses should significantly exceed 
that value. Again, GUT relations among gaugino 
masses cause the above constraint to provide a 
much more stringent (although still only indica- 
tive) bound of about 150—200 GeV. Remarkably, 
in the scenario with additional GUT unification 
of spin-zero mass parameters, just such a typi- 
cal upper bound derives from the cosmological 
constraint fl x h 2 < 1. Only in a relatively nar- 
row range of parameters where the neutralino be- 
comes nearly mass degenerate with a SUSY part- 
ner of the tau, tr, the effect of the neutralino's 
co-annihilation opens up the allowed range of m x 
up to about 600 GeV @. 

Summarizing, we can see that, despite the com- 
plexity of the neutralino parameter space and a 
large number of neutralino annihilation channels, 
one can, remarkably, select the gaugino-like neu- 
tralino in the mass range between roughly 30 and 
150 GeV as a natural and attractive DM candi- 
date. Furthermore, one is able to derive relatively 
stringent conditions on the mass range of some 
sfermions, which are consistent with our basic ex- 
pectations for where SUSY might be realized. 

3. Neutralino WIMP Detection 

3.1. Predictions 

The local halo density of our Milky Way is es- 
timated at 0.3GeV/cm 3 with a factor of two or 
so uncertainty H. For neutralino WIMPs as a 
dominant component of the halo this translates 
to about 3000 LSPs with mass m x = 100 GeV 
per cubic meter. With typical velocities in the 
range of a few hundred km/ s, the resulting flux 
of WIMPs is actually quite large 

$ = vp x /m x « 10 9 (100 GeV/m x ) x's/m 2 /s. (2) 

A massive experimental WIMP search pro- 
gramme has been developed during the last few 
years. Although a variety of techniques have been 
explored, most of them follow one of the two ba- 
sic strategies. One can look for DM neutralinos 
directly, through the halo WIMP elastic scatter- 
ing off nuclei, %A — > ^JV, in a detector. Indi- 
rect searches look for traces of decays of WIMP 
pair-annihilation products. One promising way 
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is to look for multi-GeV energy neutrinos com- 
ing from the Sun and/or the core of the Earth. 
One can also look for monochromatic photons, 
positrons or antiprotons produced in WIMP pair- 
annihilation in a Galactic halo. I will not discuss 
these indirect methods here. Perhaps I should 
only mention about an interesting new way of 
looking for WIMPs at the Galactic center p0[ . 
If there is a super-massive black hole there (for 
which there is now some evidence), it will ac- 
crete WIMPs and thus increase their density in 
the core. They will then be annihilating much 
more effectively and the resulting flux of neutri- 
nos, photons and other products from the cen- 
ter of the Milky Way may be strongly enhanced, 
even up to a factor of 10 5 in halo models with 
central spikes in their profiles. Such spiked halo 
models have been obtained in recent N-body sim- 
ulations Jig ]. 

In the following I would like to make several 
comments about direct detection and in particu- 
lar about an intriguing claim made by the DAMA 
Collaboration regarding a possible evidence for a 
WIMP signal in their data. For a more detailed 
review, see Ref. |jT9f . First, I will briefly summa- 
rize the theoretical aspects and predictions. 

In direct searches one of the most sig- 
nificant quantities is the event rate R ~ 
a(xN) (p x /m x v) (1/mjy) - the product of the 
elastic cross section a(xN) of neutralinos from 
nuclei, their flux p x /m x v and the density of tar- 
get nuclei with mass m^. 

The elastic cross section a(xN) of relic WIMPs 
scattering off nuclei in the detector depends on 
the individual cross sections of the WIMP scat- 
tering off constituent quarks and gluons. For 
non-relativistic Majorana particles, these can be 
divided into two separate types. The coherent 
part described by an effective scalar coupling be- 
tween the WIMP and the nucleus is proportional 
to the number of nucleons in the nucleus. It re- 
ceives a tree-level contribution from scattering off 
quarks, \1 ~ * X<7j as described by a Lagrangian 
C ~ (xx) (qq)- The incoherent component of the 
WIMP-nucleus cross section results from an axial 
current interaction of a WIMP with constituent 
quarks, given by £ - (X7 M 7sx) (<Z7m75<?), and 
couples the spin of the WIMP to the total spin of 



the nucleus. 

The differential cross section for a WIMP scat- 
tering off a nucleus X% with mass tua is therefore 
given by 

^^scalar ^^.axial 



da 



d\qf 



(3) 



where the transferred momentum q = mA . T " x v 
depends on the velocity v of the incident 
WIMP. The effective WIMP-nucleon cross sec- 
tions a scalar and a axial are computed by evalu- 
ating nucleonic matrix elements of corresponding 
WIMP-quark and WIMP-gluon interaction oper- 
ators. 

In the scalar part contributions from individual 
nucleons in the nucleus add coherently and the fi- 
nite size effects are accounted for by including the 
scalar nuclear form factor F(q). (The effective in- 
teraction in general also includes tensor compo- 
nents but the relevant nucleonic matrix elements 
can be expanded in the low momentum-transfer 
limit in terms of the nucleoli four-momentum 
and the quark (gluon) parton distribution func- 
tion. As a result, the non-relativistic WIMP- 
nucleon Lagrangian contains only scalar interac- 
tion terms.) The differential cross section for the 
scalar part then takes the form H 



da 



scalar 



m , -^lZf P + (A-Z)f n }*F 2 (q), (4) 

where f p and /„ are the effective WIMP couplings 
to protons and neutrons, respectively, and typi- 
cally /„ « f p . Explicit expressions for the case 
of the supersymmetric neutralino can be found, 
e.g., in the Appendix of Ref. |21j| . 

The effective axial WIMP coupling to the nu- 
cleus depends on the spin content of the nucleon 
A.q p n and the overall expectation value of the nu- 
cleon group spin in the nucleus < S p ^ n >. For 
a nucleus with a total angular momentum J we 
have 



da axv 



TTV 



2 A 2 J(J+ l)S{q), 



(5) 



with A = ■j[a p (S p )+a n (S n )). The axial couplings 
d q Aq^, a n = ^=y^d q Aq( n \6) 



a 



1 ^ 



P 



v/2 



u.d,. 



v/2 V 

u.a,; 
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are determined by the experimental values of the 
spin constants Au& = Ad^ = 0.78, Ad^ = 
AuW = -0.5 and As<» = AsW = -0.16. The 
effective couplings d q depend on the WIMP prop- 
erties and for the neutralino they can be found, 
e.g., in the Appendix of Ref. |gjj| . 

In translating a(xq) and a{xg) into the WIMP- 
nucleon cross section in Eq. (0) several uncertain- 
ties arise. The nucleonic matrix element coeffi- 
cients for the scalar interaction are not precisely 
known. Also, the spin content of the nucleon and 
the expectation values of the proton (neutron) 
group spin in a particular nucleus are fraught with 
significant uncertainty and nuclear model depen- 
dence. These ambiguities have to be considered in 
numerical calculations. Finally, in order to obtain 
(t(\N), models of nuclear wave functions must be 
used. The scalar nuclear form factor reflects the 
mass density distribution in the nucleus 

The resulting cross section for scalar (or coher- 
ent) interactions is 



2 2 

r»™^( x N) ~ G%- 

(m N + m x y 



A 



(7) 



and is proportional to the mass of the nucleus. 
For axial (or incoherent) interactions one finds 



m x m N 



{m N + m x Y 



(8) 



which can be shown to be proportional to the spin 
of the nucleus. (Gf is the Fermi constant.) 

Unfortunately, in supersymmetric models ac- 
tual calculations produce a rather broad range of 
values, R ~ 10~ 5 — 10 events/fcg/day. The rates 
are also very small. The reason why this is so is 
clear: this is because of smallness of a(xN). The 
elastic cross section is related by crossing sym- 
metry to the neutralino annihilation cross section 
cr ann which is of weak interaction strength and 
such as to give Q x h 2 ~ 1, and therefore very 
small. 

Such small event rates are clearly an enormous 
challenge to experimentalists aiming to search for 
dark matter. One may realistically expect that 
continuing SUSY searches in high energy accel- 
erators and improving measurements of Ocdm 
and the Hubble parameter will cause those broad 
ranges of R to gradually shrink. 



The not so good news is that the choices of 
SUSY parameters for which one finds the favored 
range of the relic abundance, 0.1 <Q. x h 2 < 0.15, 
correspond to rather low values of the event rate 
R < 10~ 2 events/ kg /day or so, typically about 
an order of magnitude below the reach of today's 
detectors. 

3.2. Annual Modulation 

One promising way of detecting a WIMP is to 
look for yearly time variation in the measured en- 



ergy spectrum. It has been pointed out 22 2q ] 
that a halo WIMP signal should show a peri- 
odic effect due to the Sun's motion through the 
Galactic halo, combined with the Earth's rotation 
around the Sun. The peaks of the effect are on 
the 2nd of June and half a year later. 

The effect, called "annual modulation", would 
provide a convincing halo WIMP signal. Unfortu- 
nately, in SUSY models the effect is usually small, 
AR < 5% With the absolute event rates 

being already very small, it is going to be a great 
challenge to detect the effect. 

Here I would like to make some comments 
about possible evidence for a WIMP signal in an- 
nual modulation that has been claimed by the 
DAMA Collaboration. Based on the statistics of 
14,962 kgxday of data collected in a Nal detec- 
tor over a period from November '96 to July '97 
(part of run II), the Collaboration has reported 
p4| a statistically significant effect which could 
be caused by an annual modulation signal due to 
a WIMP with mass m x and WIMP-proton cross 
section o p given as 



59 GeVl 14 GeV, 



7.0; 



-0.4 



x 10~ 6 pb 



(9) 
(10) 



at 99.6% CL, where £0.3 = Px/Po.3 stands for the 
local WIMP mass density p x normalized to pp 3 = 
0.3GeV/cm 3 . (See also Figure 6 in Ref. g]| for 
a 2(j signal region in the (m x , £o.30>) plane.) Ac- 
cording to DAMA, the new analysis is consis- 
tent with and confirms the Collaboration's earlier 



hint j 25 j for the presence of the signal based on 
4, 549 kgxday of data. 

The claimed effect comes from a few lowest 
bins of the scintillation light energy, just above 
the software threshold of 2 keV, and predomi- 
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nantly from the first bin (2 — 3 keV). This is 
indeed what in principle one should expect from 
the annual modulation effect. DAM A appears 
confident about the presence of the effect in their 
data, and claims to have ruled out other possi- 
ble explanations, like temperature effects, radon 
contamination or nitrogen impurities. Accord- 
ing to DAMA, the effect is caused by single hit 
events (characteristic of WIMPs unlike neutron 
or gamma background) with proper modulation 
of about one year, peak around June, and small 
enough amplitude of the time dependent part of 
the signal. 

Nevertheless, several experimental questions 
remain and cast much doubt on the validity of 
the claim. Here I will quote some which I find 
particularly important to clarify. First, as stated 
above, the claimed effect comes from the low- 
est one or two energy bins. This is indeed what 
one should expect from an annual modulation sig- 
nal. But is the effect caused by just one or two 
energy bins statistically significant? This is es- 
pecially important in light of the fact that the 
shape of the differential energy spectrum dR/dE 
in the crucial lowest energy bins as measured 
by DAMA is very different from the one mea- 
sured by Gerbier, et al. |2q| fo r the same detec- 
tor material (Nal). In Ref. []26f the corrected-for- 
efficiency dR/dE is about 10 events/ kg /keV/ day 
at 3 keV, decreasing monotonically down to about 
2 events / kg /keV/ day at 6keV. (See Fig. 15 
in Ref. Q.) In contrast, DAMA's spectrum 
shows a dip down to 1 events/fc<?/fceVyday 
at 2keV, above which it increases to nearly 
2 events / kg /keV/ day at 4keV @|7)- :t is abso- 
lutely essential for the controversy of the shape of 
dR/dE in the lowest energy bins to be resolved. 
Furthermore, DAMA's data from the second run 
(~ 15000 kg x day) shows that the background in 
the crucial lowest bin ([2,3] keV) is only about 
half or less of that in the next bins [^7j . One may 
wonder why this would be the case. Examining 
more closely the data in the constituent nine Nal 
crystals, one finds a rather big spread in the event 
rates [g7J. I n detector 8, in the lowest energy bin 
one finds no contribution from the background 
whatsoever! 

Two other groups which have also used Nal 



have reported [|26|,f28| robust evidence of events of 
unexpected characteristics and unknown origin. 
The data of both teams has been analyzed using 
a pulse shape analysis (PSA). A small but statis- 
tically significant component was found with the 
decay time even shorter than the one expected 
from WIMPs or neutrons. While the popula- 
tion of those events appears to be too small to 
explain DAMA's effect, a question remains not 
only about their origin (contamination?, external 
effect?) but also how they contribute to the en- 
ergy spectrum in the crucial lowest bins. DAMA 
claims not to have seen such events. 

DAMA has not yet published the data over 
a full annual cycle|] In particular, no evidence 
has been published of the signal going down with 
time. So far, the claimed signal has been based on 
two periods of data taking. Ref. |2f| was based 
on two relatively short runs, one in winter and 
one in summer. The second analysis |24| ] used 
the data collected between November and July. 
Much more data has been collected and will soon 
be published. One should hope that a full and 
clear analysis of statistical and systematic errors 
will also be performed. 

In Ref. [p9[ annual modulation was reanalyzed 
for germanium and Nal. It was concluded that the 
effect would be too small to be seen with current 
sensitivity. Particularly illuminating is Fig. 6. a 
where DAMA's data from Ref. |25| (run one) was 
re-plotted along with an expected signal for the 
modulated part of the spectrum for the central 
values of the ranges of the WIMP mass and cross 
section (a p ) selected by DAMA. One can hardly 
see any correlation between the data and the ex- 
pected signal. 

A controversy over DAMA's alleged signal is 
of experimental nature. One may therefore hope 
that it will be definitively resolved. It would be 
of particular importance for another experiment 
using different detector material to put DAMA's 
claim to test. The CDMS cryogenic detection ex- 
periment using germanium and silicon crystals at 
Stanford has now reached an adequate sensitivity 
and has already ruled out about a half of the re- 
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WIMP Mass [GeV] 



Figure 1. The current and some recent limits on 
the WIMP-proton spin-independent (scalar) cross- 
section. The legend is as follows: dot: Heidelberg- 
Moscow Ge @, solid: DAMA Nal |l|, dash: CDMS 
1999 4 kg x day Ge with neutron background sub- 
traction Q . Marked also are the 2 - a DAMA re- 
gions: solid closed curve: 15,000 kg x day, filled re- 
gion: 20,000 kgxday E3|. 



gion selected by DAMA (27) .| Let us hope that 
DAMA's claim will be falsified soon. 

Figure [l] shows the current status of WIMP 
searches in the plane of the WIMP-proton scalar 
cross-section o v versus the WIMP mass.^j A 
claimed DAMA region is also indicated. 

It is worth pointing out that the claimed sig- 
nal region selected by DAMA was too re- 
strictive as it did not include astrophysical un- 
certainties. The effect is rather sensitive to as- 
sumptions about a model of the Galactic halo. 
In the DAMA analysis the peak of the WIMP 



4 See Note Added. 

5 The figure has been produced with the help of 
a particularly useful plotting facility of Gaitskcll 



(Maxwellian) velocity distribution was fixed at 
i>o = 220 km/ s. Since the Galactic halo has not 
been directly observed and there still is a consid- 
erable disagreement about a correct halo model, 
quoted error bars for vq and the local halo density 
should, in my opinion, be treated with much cau- 
tion. Varying vq within a reasonable range leads 
to a significant enlargement of the region selected 
by DAMA, H as first shown in Rcf. §|. This 
is so because of the significant dependence on vo 
of the differential event rate spectrum, as can be 



Mandic which is now availa 



http: / /cdms. berkeley.edu/limitplots/ 



ble on the Web: 



seen in Fig. 1 in Ref. |32|. 

As a result, the upper limit of m x correspond- 
ing to the claimed signal region increases consid- 
erably, from about 100 GeV for the initially as- 
sumed value «o = 220 km/ s up to over 180 GeV 
at v = 170 km/s Q (at 2er CL). On the other 
hand the range of £o.30> is not much affected. 

Assuming that the effect claimed by DAMA 
were indeed caused by DM WIMPs, it is interest- 
ing to ask what ranges of SUSY parameters and 
LSP relic abundance it would correspond to. This 
issue was addressed by three groups, each using 
their own code for computing the direct detec- 
tion rate and the relic abundance |34|-|3(|. While 
there is some difference in approach (somewhat 
different ranges of input parameters, techniques 
and assumptions, etc.) the overall outcome is 
that it is indeed possible to find such SUSY con- 
figurations which could reproduce the signal but 
only for large enough tan/? (typically above 10 
although SUSY points with smaller values can 
also sometimes be found). More importantly, the 
corresponding values of Sh,h 2 are rather small, 
typically below 0.06 [p4|-p6[ although somewhat 
larger values can also sometimes be found. These 
are clearly small values, well below the favored 
range 0.1 < fl x h 2 < 0.15. This is illustrated in 
Figure | f§. 

However, one should remember that these re- 
sults have been obtained using commonly used 
but often rather uncertain values of several input 
parameters. In addition to astrophysical ones, the 
nuclear physics and quark mass inputs in calcu- 
lating the scalar cross-section for the neutralino- 
nucleus elastic scattering are rather poorly de- 
termined, as mentioned above. The effect of the 
latter has been recently re-analyzed in Ref. B7]] 
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H<0, tan/S= 30, v = 220 km/s, 0.3<f 03 <3 
0.0001 p — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r- 




WIMP Mass [GeV] 



Figure 2. A scan of SUSY points in the plane 
(m x ,£o.30>) where £0.3 = p x /(0.3 GeV /cm 3 ). The 
legend: circles: 0.1 < Q x h 2 < 0.15 (cosmologically fa- 
vored range), crosses: 0, x h 2 > 0.25 (excluded), small 
points: 0.02 < Q x h 2 < 0.25 (conservative range). 
The 2 — o region (upper left) of DAMA is marked 
with full thick dots. 



and found to affect the overall scalar cross-section 
by a factor of ten or so. 

It is worth presenting Figure |3| again with an 
outlook for some expected limits (as claimed by 
the respective groups) and compare them with 
predictions of minimal SUSY obtained for very 
broad ranges of SUSY parameters and neglect- 
ing nuclear input uncertainties mentioned above. 
The SUSY region is bounded from above by a 
somewhat arbitrary requirement fl x h 2 > 0.025. 
From below it is limited by a generous bound 
Vt x h 2 < 1. The SUSY points falling into the cur- 
rently expected range 0.1 < il x h 2 < 0.15 (not 
indicated in the plane) form a sub-region reach- 
ing up to roughly a few times 10 -6 pb. 

It is clear that today's experiments are now 
only reaching the sensitivity required to begin 



Figure 3. The current and some future reach lim- 
its on the WIMP-proton spin-independent (scalar) 
cross-section compared with conservative predictions 
from minimal SUSY (large shaded region). The leg- 
end is as follows: upper dot: Heidelberg-Moscow 
Ge Q, solid: DAMA Nal Q, upper dash: CDMS 
1999 4 kg x day Ge with neutron background subtrac- 
tion J27]] . Marked also are the 2 — a DAMA regions: 
solid: 15,000 kgxday, filled: 20,000 kgxday §4|. 
Future reach of some experiments (as claimed by 
the respective groups): lower dash: CDMS at 
Soudan, lower dot: CRESST. Some other experi- 
ments (e.g., Edelweiss, UKDMC-Xe, GENIUS) ex- 
pect to reach roughly similar sensiti vities. (Source: 
[http:/ /cdms.berkeley.edu/limitplots/). 



testing predictions coming from SUSY. What I 
find promising is that several experiments us- 
ing different detector materials and often differ- 
ent methods of [attempts at] distinguishing signal 
from background will explore a large fraction of 
the SUSY parameter space within the next few 
years. Especially reassuring would be an observa- 
tion of a positive signal in more than type of DM 
detector, although many experimentalists would 
probably remark that I am asking for too much. 
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Time will tell. 
4. Axinos 

Axinos are a natural prediction of the Peccei- 
Quinn solution to the strong CP-problem and 
SUSY. The axino is the fermionic partner of the 
axion. Similarly to the axion, the axino couples 
to ordinary matter with a very tiny coupling pro- 
portional to l/f a where 10 9 - 10 10 GeV < f a < 
10 12 GeV. 

It is plausible to consider the axino as the LSP 
since its mass is basically a free parameter which 
can only be determined in specific models. As we 
have seen above, the neutralino has been accepted 
in the literature as a "canonical" candidate for 
the LSP and dark matter. But with current LEP 
bounds between about 30 and 60 GeV (depend- 
ing on a SUSY model), it becomes increasingly 
plausible that there may well be another SUSY 
particle which will be lighter than the neutralino, 
and therefore a candidate for the LSP and dark 
matter. 

Primordial axinos decouple from the thermal 
soup very early, around T ~ f a , similarly to the 
axions. The early study of Ragagopal, et al. |Q 
concluded that, in order to satisfy Vlh 2 < 1, the 
primordial axinos had to be light (< 2 keV), cor- 
responding to warm dark matter, unless inflation 
would be invoked to dilute their abundance. In 
either case, one did not end up with axino as cold 
DM. 

However, it has recently been shown |^] that 
the axino can be a plausible cold dark matter can- 
didate, and that its relic density can naturally be 
of order the critical density. The axino can be 
produced as a non-thermal relic in the decays of 
heavier SUSY particles. Because its coupling is 
so much weaker, superparticles first cascade de- 
cay to the next lightest SUSY partner (NLSP) 
for which the most natural candidate would be 
the neutralino. The neutralino then freezes out 
from thermal equilibrium at Tf ~ m x /20. If it 
were the LSP, its co-moving number density af- 
ter freeze-out would remain nearly constant. In 
the scenario of Covi et al. (CKR) ||, the neu- 
tralino, after decoupling from the thermal equi- 
librium, subsequently decays into the axino via, 



e.g., the process 

X^aj (11) 

as shown in Fig. 1 in Ref. ||. This process was 
already considered early on in Ref. |39| (see also 
Ref. J3^]) in the limit of a photino NLSP and only 
for both the photino and axino masses assumed 
to be very low, m~ < 1 GeV and m~ < 300 eV, 
the former case now excluded by experiment. In 
that case, the photino lifetime was typically much 
larger than 1 second thus normally causing de- 
struction of primordial deuterium from Big Bang 
nucleosynthesis (BBN) by the energetic photon. 
Avoiding this led to lower bounds on the mass 
of the photino, as a function of f a , in the MeV 
range Q. 

Because both the NLSP neutralino and the 
CKR axino are both massive (GeV mass range), 
the decay ([ll]) is now typically very fast. In the 
theoretically most favored case of a nearly pure 
bino, |ll|Jl5tl the neutralino lifetime can be writ- 
ten as 

r " 3 - 3xl ° b { io" Gev ) \r^r) (12) 

where the factor NC a yy is of order one. One 
can see that it is not difficult to ensure that 
the decay takes place well before 1 second in 
order for avoid problems with destroying suc- 
cessful predictions of Big Bang nucleosynthesis. 
The axino number density is equal to that of the 
NLSP neutralino. Therefore its relic abundance 
is n~h 2 = (m~/m x ) U x h 2 . The axinos are ini- 
tially relativistic but, by the time of matter dom- 
inance they become red-shifted by the expansion 
and become cold DM. 

There are other possible production mecha- 
nisms and cosmological scenarios for massive ax- 
inos J^j Even if the primordial population of axi- 
nos is inflated away (which would happen if the 
reheating temperature T ro h <C f a ), they can be 
regenerated from thermal background processes 
at high enough T re h. 



For a recent comprehensive study, see 
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5. Conclusions 

Looking for the invisible is not easy but is cer- 
tainly worthwhile. The discovery of dark matter 
will not only resolve the mystery of its nature but 
is also likely to provide us with much information 
about the particle physics beyond the Standard 
Model. What I find most encouraging is that we 
have a very good chance of fully testing the neu- 
tralino as a WIMP by the end of the decade. De- 
tecting other candidates, like the axino, will be 
the task for the more distant future. 



Note Added: 

After the Conference both the DAMA 
and CDMS Collaborations published new re- 
sults. Based on the combined statistics of 
57, 986 kgxday of data collected in a Nal detec- 
tor since November '96, the DAMA Collaboration 
reported El[ a statistically significant (4cr CL) 
effect which it interprets as being caused by a 
WIMP annual modulation signal. 

The CDMS experiment using germanium and 
silicon crystals at Stanford published jEl a new 
limit on scalar WIMP-proton cross section. The 
new result is based on the total of 10.6 kgxday of 
data collected a current shallow site (17 mwe) at 
Stanford during 1999. A powerful event-by-event 
discrimination method allows CDMS to reach a 
sensitivity matching that of DAMA with only less 
than 0.2% of DAMA's statistics. 

The new 90% CL CDMS limit excludes most 
of the signal region claimed by DAMA. In par- 
ticular, it fully rules out the previous 2 a region 
based on the combined data of 19,511 kgxday 
from runs I and II [^ij and the new 3 a region at 
more than 84% CL. 

An updated compilation of these and 
other data can be found on the Web: 



3. L. Covi, L. Roszkowski and J.E. Kim, Phys. 



bttp: / /cdms. berkeley.edu/limitplots / . ) 
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